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^ ' Abstract 

■ The chiral symmetry breaking in the 4-dimensional QED with the chirally 

! invariant four-fermion interaction is discussed by using a novel path integral ex- 

pression in terms of the field-strength tensor. In the local potential approxima- 
tion, we find that the chiral symmetry is spontaneously broken for any nonzero 
Oh! gauge and four-fermion couplings on the tree level of an auxiliary field a. The 

(— I I present approach allows us to easily include higher orders of the gauge coupling 

so that the effective potential up to the sixth order is obtained. 

• 1— I 

X. 

c5 ■ 1 Introduction 

The standard model of elementary particles has had good agreement with experiments, 
in which masses of quarks and leptons are generated by the Yukawa interactions with 
the Higgs bosons. There, however, have been no experimental observations of the Higgs 
boson, and moreover, there is a theoretical problem of the "naturalness." Therefore 
models without any elementary scalars have been considered [||, |^ as effective field 
theories to understand the origin of fermion masses: Those are generated dynamically 
through the interaction 

^ + (?^75^)2] , (1.1) 

which is suppressed by inverse powers of a high energy scale A, at energies much below 
it, that is, g'^A'^ ~ 0{1). It is called the Nambu-Jona-Lasinio (NJL) model[Q and has 
been generalized to the one with the gauge interaction |^: 

C = -\f^.F^^ + iii^^id^ - ieA^)ij + ^ [i^ijf + (^^75^)'] • (1.2) 
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Although, in this model, nonperturbative approaches such as the Schwinger-Dyson 
equation with the ladder approximation!^, |^ or renormalization group give a non- 
trivial phase structure for chiral symmetry in a plane of gauge and four-fermi cou- 
plings, they suffer from gauge dependence. In the perturbative framework the inversion 
method^ does give a gauge-independent result ||], higher order calculations are tough 
tasks. 

If we employ the perturbative effective action approach to the problem of chiral 
symmetry breaking(xSB), there is some room for improvement or simplification of the 
calculation from two viewpoints: (i) In view of the nonperturbative renormalization 
group approach]^, the critical behavior is successfully described in the local potential 
approximation: taking the lowest order of derivative expansion in the effective action. 
We, therefore, expect that it could also be good in the perturbative calculation, (ii) 
After integrating with respect to matter (fermionic) fields the result must be gauge 
invariant, that is, be written as a functional of the field-strength tensor F^^. Therefore 
it is preferable to rewrite the functional integral of the gauge potential into that of 
F^y. This could be a great help to the gauge-invariant calculation. 

Indeed in the lattice gauge theory, a change of variables from gauge potentials to 
field strengths has been given[0], whose recipe, however, is specialized to the lattice. 
In the continuum theory, the field-strength formulation has already been proposed ||TT . 
It is, however, based on the special choice of gauge — the coordinate gauge x^A^^ = 0, 
which cannot fix the gauge at the origin and moreover yields complicated fermionic 
currents with line integrals in the lagrangian. 

In this article we give a perturbative and gauge-independent method for calculat- 
ing the dynamical fermion mass under the local potential approximation. Inclusion 
of higher orders of the gauge coupling constant is simpler in this method. We use 
the functional integration not as a method simply reproducing perturbative diagrams 
systematically, but as "integration" with nice features: Changing variables and intro- 
ducing auxiliary variables can easily be made. We first construct a Euclidean path 
integral expression in terms of the gauge field strength for the gauge sector with a 
conserved source J^: We start from the canonical formalism for quantization and find 
a suitable change of variables, which is the contents of Sec. |[ Second, in Sec. ^ the 
fermionic partition function, minimally coupled to the gauge potential, in the local 
potential approximation is obtained by introducing auxiliary fields and by utilizing the 
Fock-Schwinger proper-time method [|T^, At last combining the above two partition 
functions and integrating with respect to the gauge field strength to obtain an effective 
potential of the dynamical fermion mass. The final Sec. ^ is devoted to a discussion. 



2 Functional Integral in terms of the Field Strength 

In this section we construct a Euclidean path integral expression in 4-dimensional 
abelian gauge theory coupled to a conserved current J^, based on the canonical for- 



malism 14 



C = -^F,,F'^'' + J^A>^. (2.1) 
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In order to show the gauge independence manifestly we fix the gauge in terms of 
an arbitrary real function satisfying fl^, |15| 



the gauge fixing condition is|14 



d,r{x) = 5\x)- (2.2) 



^d''y<P,{x-y)A^{y) = Q, (2.3) 
which is satisfied by 

A^^{x-(t>) = A^{x)-d^^ [dSUx-y)A^{y), (2.4) 



for an arbitrary A^^{x). Here we take a (f)fj,{x) whose support is 3-dimensional and 
spacelike: 

r{x) = (0,f (x)5(xo)), rip) = (0,r(p)) (2.5) 

with 

Vj\x) = S^x), p9{p) = -^■ (2.6) 



If it were not for this restriction, A^{x; 0) in Eq. ( |2.4| ) is nonlocal in time, so that we 
cannot follow the canonical procedure. (Here and hereafter repeated indices i, j, k, 
etc., imply the summation over 1 to 3.) The relation ( ^.41 ) turns into 

A%x;(f)) = A%x) + f{-iV)A^{x), (2.7) 

A\x;(f)) = A\x) -V'fi-tV)A^{x), (2.8) 

where 

f{-zV)A^{x) ^ l-0-^e^P-y^{p)A^{p) = jd'yfix - y)A\y) (2.9) 

and the time argument has been omitted. Hereafter we use an abbreviation like 
Eq. (|2.9| ) for simplicity; all the functions of — iV including 1/|V| and 1/V should 
be understood in the momentum space representation. 

The gauge-fixed variable A^{x](j)) must have two degrees of freedom out of three 
y4*(a;)'s, which can be singled out by considering the norm of the functional space 

jd^x5A\x-(t))5A\x-(l)) = jd^x5A^ix)M^''{-iV)5A\x), (2.10) 

where 

M^'ip) ^ + *(p)/ - iTPf\p) + p'f*{p)f\p). (2.11) 
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This matrix can be diagonalized as 



( 1 



where n^^^^'s are given by 



a/3 



1,^{p) ^ tp' + p'f{p) Np\p'\f{p)?-i) 



n%){p) = Inl- 
and form an orthonormal base obeying 
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Q = l 



In view of Eqs. ( p.lO| ) and ( p.l2| ), genuine physical variables are given as 
A^-\x) ^ nl^{-^V)A'{x) = nf,)(-zV)A'=(a;; 0), (a = 1,2). 
With these variables the action is 

S = Jd^'x !^~[d^A,{x-(j))-d,A^{x;^)f + J^{x)A>^{x-^) 



d\ \IJ2a^''\x)V'A^''\x) + i \A^'\ 



1 2 



a=l 



- \M'\x)V'\f{-iV)\^A^'\x) + \ [VA^x-A)]- 



jQ{x)A^{x-(t)) - J{x)-[A{x](t) 



where 0)] in the last line is given by 

A\x-<P) 



n;t)(-^V)A«(a;) 



+ 



nf2) (-^ V) + i-zV) V - | /(-^ V) P - 1 



'-(3) 



^(2)( 



The Hamiltonian is written as 



H{t) = /rf'a; j - ^ ( [n(") (x)] ' + [Vv4(") (a; 



;v^(v^l/Hv)l^ + i) 



+ 2'^o(a:)|/(-zV)|Vo(x) + J(x)-[A(a;; 



in terms of four dynamical variables, two A^°'\x)^s and their canonical conjugate mo- 
menta 



n«(a;) 



v^(v^|/HV)|2 + i) 



Joit,x), 



(2.19) 
(2.20) 



where A^{x; </>) has been eliminated by 



V^(V^|/(-^V)|2 + l) 



The field strengths are given as 



U^'\x) + 



1 



l/(-^V) 



-A°(a3;0) + Jo(t,a3) = 0. (2.21) 



E' = F 



iO 



-n 



(1)' 



v)n« - nf2)(-^v)n(2) + r (-zv) Jo, 



1 



jk 



|V|n5i)(-iV)A(2) - |V|nj2)(-«V)A( 



1) 



(2.22) 



(2.23) 



Following the standard canonical procedure we obtain the path integral representation 
of the partition function Zt[J] = Tr(e~"^^)inT 



Zt[J] 



exp 

-5t([n' 



r, X 



! 'fxAi n(°) (r, a;)i(") (r, a;) 



r, a; 



+z J4(r, £c) 2 n^^^ (r, a;) 



+ -J4(r, a;)|/(-zV)|V4(r, a;) - J(r, a;)-[A(r, x; 



(2.24) 



where [A(r, a;; 0)] is given by Eq. ( p.l7|) with A{x] (p) A(r, x; 



d X, 



dr / d X, Ja = iJn 



(2.25) 



and the periodic boundary condition yl'-"''(T, a;; 0) = A*^"''(0, a;; 0) should be under- 
stood. 

Our purpose is to rewrite the expression (|2.24]) of four variables 11" and A" into 
that of six variables and B^. In view of Eq. ( p.l4|) quantities which are proportional 
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tV) are missing in their expressions ( p.22| ) and ( p.23| ). To this 
end, we introduce e and j3 with the aid of the delta functions in the functional measure, 



to n^^g^f— iV) or n 



(3)^ 



E' = 
B' = 
or equivalently 



fDHV)n(i) _ ^» (_,v)n(^^ + n^^%{-tV}e - tr{-tV)J,, (2.26) 



r(2) ^ n^* 



|V|4)(-zV)A(2) - |v|n^(2)HV)A« + ni3)HV)/3, 



n 



(2.27) 

(2.28) 
(2.29) 



d^x5E'5E' = yd='x((5n(^)(5n(^) + 5n(2)^n(2) + fefe), (2.30) 

jd^x5B'SB' = jd^xi-5A^^^V^SA^^^ -SA^^^V^SA^^^ +5P5P), (2.31) 



P = nX,){-iV)B^ 



The norm of the functional space is given by 



so that 



25n(")pA(")p£D/3 Jj5(e)5(/3) = VE'VB'\\_5{y'E' - iJ^)5{y'B'). 

X X 

We then arrived at the desired result 
Zt[J] = fvE'VB'YlSiV'E' -iJ4)5iV'B') 

X 

X exp jd^x^!^iE\T,x)e'^''^B''{T,x) - ^ (^[E'\T,x)f + [B' 



(2.32) 



(r, X 



+ J\T,x)e^^'^^B>^iT,x)^ 



(2.33) 



which is apparently free from the choice of f^{x) as expected. 

To carry out the integration in the next section, it is convenient to introduce E^ as 



so that 

giving 
Zt[J] 



E = E^ — i — ttJa, 



(2.34) 
(2.35) 



VE'VB' JJ 5{V'E")5{V'B') exp 



[E{t,x)\ +[B\t,x 
+ J\t, x)e'^''^B\T,x)'^ 



d*x^\iE^{T,x)e'^''^B''{T,x) 



1 



-J4{t,x)^J4{t, x) 



(2.36) 
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Further introducing an auxiliary field p\ 

1 



Vp exp 



so as to cancel the ( ^4)^ term in Eq. ( p2.36| ), decomposing E and as Q 



(2.37) 



(2.38) 



Q=l a=l 

and integrating with respect to e^'s and /?3, we obtain 
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Zr[J] = /p/5„r'p[Det|V|]^2exp 



-ip(r, x)^ - J4(r, x)-^p(r, a;) 



4e 



H 2 



r/3Q(r,a;) 



+ /5„(r,a;) 



(2.39) 



This form with further changes of variables is used in the next section, where we regard 
the coefficients of the sources in Eq. (p.39|) as gauge potentials: 



A, = 
A, = 



1 



1 



[ng)(-^V)/3i-nJ*)(-zV)/32]. 



(2.40) 
(2.41) 



Therefore field strengths are given as 



i4 



Vi - Ai 



Fij = 'ViAj — VjAj 



eijfc [nfi)(-iV)A +nf2)(-2V)/32] 



so that 



.a=l 



(2.42) 
(2.43) 

(2.44) 



^Any choice for the complete orthonormal base {rt^^^(— iV)} can be taken. Here, we only assume 
the same nj3j(— iV) as Eq. ( 2.13 ). 
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3 X^SB in QED4 with the chirally invariant four- 
fermion interaction 



In this section, we consider fermionic system coupled minimally to the "gauge po- 



tentials" (|2.40|) and (|2.41|) with the chirally invariant four-fermion interaction. The 
partition function is 



Z[A] = VipViljexp 



(3.1) 



where the gauge coupling constant has been absorbed into A^. Our scenario is as 
follows: introduce auxiliary fields to cancel the four-fermion interaction, and then 
integrate with respect to the fermionic fields and finally the gauge field strengths with 
the aid of the representation in the previous section. The result is a tree potential of 
the auxiliary field, with which we examine the dynamical mass generation of fermions. 

After introducing auxiliary fields a and vr, as usual, and integrating with respect to 
if) and ?/', we have 



Z[A] = VaVnexp 



— Id 



+ InDet [7^(9^ - iA^}, + (x + i^^-n] 



(3.2) 



Shifting as cr — > m + 0"' and vr ^ vr' and ignoring a' and tt', we obtain the tree level 
potential of oim): 



Z[A]q = exp 



d'^x^ ^ + In Det - lA^) + m] 



(3.3) 



The exponent of Eq. (|3.3| ) must be a functional of the field strength F rather than 
as a regularization preserves the gauge invariance. 

We here employ the local potential approximation: We adopt the lowest order of 
derivative expansion, that is, discard any terms with differentials like F^JUj^Ffj^^y, to ob- 
tain a polynomial of F^^/. This approximation seems to be valid, since we are interested 
only in a low energy phenomena, the chiral symmetry breaking, where contributions 
from large should be much less important. 

The functional form of the effective action under this approximation can be obtained 
nonperturbatively by the Fock-Schwinger's proper time method|]T2|, pr5| : 



Z[A]o = exp 



d^x. 



m 



+ 



2^' 2(27r) 



-lim 

— s->0 



drr' 



D_ 

' 2 ^ 



G{tF) 



(3.4) 



where 



G{F) 



F+F^ coth(F+) coth(F_ 



F F + F F 




± 



F F 



F F 



(3.5) 
(3.6) 
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To evaluate the r-integration in Eq. 

G{F) : 



we expand G{F) as 



1 + \{Fl + F')-^ [{Fl + F'f - 7FlF'] 



1 



945 



[2{Fl + F'^f - 13FlFl{Fl + F^)] + O(F^). 



(3.7) 



We need some regularization: In order to reproduce the NJL result in the limit e — 0, 
we need an ultraviolet cutoff A which is introduced by a modification of the range 
of the r-integration to [l/A^,oo). However, this proper-time cutoff breaks the gauge 
invariance in the same way as the momentum-space cutoff, contrary to the dimensional 
regularization {D = 4 — 2e). To overcome this difficulty, we employ both regularizations 
at the same time. We use the cutoff for the zeroth order in the gauge coupling, since 
it has nothing to do with gauge fields. While, for higher orders, the dimensional 
regularization is usedQ. The result is 



Z[A], 



exp 



m 
2? 



+ 



A' 



1/1 , yU^ 

+ - 3 + ln- 
3 Ve m 



F F 



F F 



+ 



45m^ 



F F 



+ 



315m8 



F F 



- 13 



F F 



2p p 



+ 0(F^ 



where 



Ei{z) 

1 

e 



-zt 



dt- 



7 + ln27r, 

e 



(3.8) 

(3.9) 
(3.10) 



and /i is a renormalization scale. [Note that Ei{z) > for any real z{> 0).] Recall 
that the gauge action is written as (— l/4e 
such that 



larelFfiuFfj.u to define a renormahzed charge 



1 



+ 



1 1 



el (/x) el (/i) 127r 



2 



(3.11) 



^If we adopt the dimensional regularization from the zeroth order we have a different critical 
coupling due to the lack of a quadratic term of the renormalization scale ^ in the effective action. 



Compare two expressions of the zeroth order effective action, Eq. ( |3.13 ) by the proper-time cutoff 
after expanding in terms of to^/A^, 



2g2 



m 



and one by the dimensional regularization, 



Vo(m) 



252 



167r2 



A4 



A2 



7 + - + In ■ 



A2 



In 27r — 7 + 



In- 
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where the first term of the right-hand side is the bare part and Z3 is the wave function 
renormahzation constant. 

Now we turn our attention to the functional integration of the gauge field strength. 
The total partition function is given by combining Eq. (|3.8|) with the result in the 
preceding section: 



Zlm] 



where 



exp 
X exp 



d X, 



1 



VPJ)p[Y)ei\V\] 



(3.12) 



1 n 
- in — - 

3 



2 ^ p P 



+ higher orders 



29' 



(3.13) 



"higher orders" are those of Eq. ( |3.8| ) with F^u F^y , and F^y are given by Eqs. ( |2.42|) 

and (HD. All variables in Eq. ( gl2D are F ni^Fnu and FnyF so we can change the 
integration variables from (3a and p to S, T, and U: 

^ ^ 2(52 + r2 + f/2), (3.14) 

2{S'^-T'^), (3.15) 



yielding 



V(3J)p[J}ei\V\ 



(3.16) 



from Eq. ( |2.44D . As for F^j^i^F^^, though it cannot be expressed by the total divergence 
any more, its expectation value in the space-time integral would vanish, 

d^x^(F^,f^,) = jd^xMS^ - T^)) = 0, (3.17) 

due to a symmetry S ^ T oi the effective action ( p.l2| ); the dependence on S and T 
appears only in the forms S"^ + and {S'^ — T^)^. 

Note that Z[m\ ( |3.12|) is a trivial product of integrals at each space-time point. We 
evaluate these integrals with the help of the WKB approximation. For a technical rea- 
son we discard 0{F^) terms in Eq. (|3.8|) . (This enables us to obtain the stationary point 
analytically. See below.) Discretizing space-time with (lattice spacing)^ = 327r^/A^ 
and neglecting the irrelevant factor Det[— □ ]~^, we obtain 



Z[x] 



n 8(32 
{site} 



3 

7r)2 



ds dt du 



[0,oo)3 



X exp 
1 



Vq[X) 



( 167r^ 4 , p? \ , •) 9 9s 



e|(/i) 3 xA2 
^3^,[4(.^ + t^ + .^)^-7(.^-tT] 
2 



315x4 



[8(^2 + t2 + ^2)3 _ 13(^2 ^ ^2 ^ ^2)(^2 _ ^2)2] | 



(3.18) 
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where the first coefficients are normahzation factors to cancel the Gaussian integration; 
dimensionless parameters, x — /A?, s — S/A^, etc. have been introduced; and 



vo[x) = = + 2(1 - x)e-^ + 2x^Ei{x). 



A4 



^2A2 



(3.19) 



Introducing a polar coordinate (r, 9, ip) with = s'^+t'^+u'^ and sin^ ^ cos 2ip = s^—t^ 
and integrating with respect to 6 and </?, we find 

/■oo 

Z[x]^ n 47r(327r)i / cir exp [- {^;o(a;) + ^;^(r; x)}] , (3.20) 



where 



{site} 



Vp{r-^x) = A{x)r^ — B{x)r^ + C{x)r^ — \nr'^ , 



IGtt^ 4 , 

- Inx, 



B{x) 
C{x) 



e|(//) 3 xA2 e2(A) 3 
32 

136 



4725^4 ' 



(3.21) 
(3.22) 

(3.23) 
(3.24) 



and O(r^) terms in the exponent has been neglected. 

Let us first consider the lowest correction, up to the O(r^) term. The r- integration 
is analytically performed to give 



Z[x] = JJ 2^5/^3 exp[-^;(x) 
{site} 



(3.25) 



where 



v{x) — Vq{x) ^ — \yiA{x) 
2 

At '\ 
= -^ + 2(l-x)e-^ + 2x2£;i(x) + Jin 
Ct 2 



47r 



4, 

- mx 
3 



(3.26) 



with G = q'^A'^/Att'^ and a{A) = e^(A)/47r. The argument of the logarithm becomes 
negative at the Landau pole x = exp[37r/Q;(A)], but we do not care such a heavy fermion 
and thus assume that the argument is always positive. The stationary point x* defined 

by 



1 dv{x) 



4 dx 

should satisfy the stability condition 

1 d'^v{x 



1 -X* *p f *N 3q;(A) 
--e J-8(3^_^(^)ina;*)a;* 







(3.27) 



4 



97ra(A)-3a(A)^(lna;* + l) 



8(37r-Q;(A)lnx*)2x 



(3.28) 
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for which x* < e ^ ~ 0.368 is a sufficient condition. To solve Eq. ( p.27| ) we set a 
condition: 

1 > -xlnx > (1 - 7)x > (3.29) 

which is fulffiled if x < 1. x 10~^. (In the actual situation m ~ 1 MeV and we 
should take A > 1 TeV, a lower bound of compositeness from experiments , to give 
X < 10~^^.) Under this condition Eq. ( |3.27|) becomes 



— — 1 — X In X 

Gr 



3a(A) 



^(37r — a{A) lnx)x 
where we have used the expansion of -E'i(x) for x <^ 1, 

Ei{x) = —7 — Inx + 0(x). 



(3.30) 



(3.31) 



For a (A) = there exists a nonvanishing solution only if G > 1; therefore the critical 
coupling Gc is 1. The solution for G ~ 1 is 



G-1 



X 



ln[G- 1]' 

For a (A) > the solution is obtained in two separate regions: (i) x ^ exp 

3G 



(3.32) 



X 



"81n(f) 



a(A)_ 

(3.33) 



which is independent of a (A), (ii) exp 



X 



a(A) 

a{A)G 
87r(l-G)' 



< X < a(A)/4; 



(3.34) 



Thus there always exists a nonvanishing solution for a given G, that is, Gc = 0. 
[Solutions at several a(A)'s are depicted in Fig. |l|.] The actual situation, a{A) ^ 1/137, 



X ~ 10 and A ~ 1 TeV, lies in the case (ii) (exp 
Eq. (KM) 



3n 



a(A) 



1.7 X 10 



-561^ 



ttG = — — ~ 1.1 X 10"*; 

Arc a (A) 



and from 



(3.35) 



g is highly suppressed even at this scale (A ~ 1 TeV). 

Now we include the higher order terms. We evaluate the r-integration with the 
WKB approximation; 



/>oo ^ 

/ dr exp [—Vp{r; x)] ~ (27r)2 exp 
Jo 



- { Vp{ro{x); x) + i In [t;^(ro(x); x)] 



(3.36) 
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Figure 1: Squared dynamical mass of fermion m? / shown as a function of the four- 
fermion couphng constant (yf^A^/47r^ = G for several fixed gauge coupling constants 
a (A) = e^ (A)/47r, obtained from Eq. ([OtD . 



where the superscript ' denotes the r-differentiation and ro(x) is a solution of the 
stationary point equation 



v'^{r- x) = - [A{xy - 2B{x)r^ + ^C^x^ - l] = 0, 



which is the cubic equation of r^. There exists only one real-positive solution for r 



(3.37) 

2 



r^^x) 



56 
153 



+ - 



3 /175 



2 V 102 



X3 [P{x) + g(x)]3 - [P{x) - Q{x)Y , 



where 
P{x) 
Q{x) 



1 + 



802816 
47403225^ 



x"^ — 



112 
153 

2 



x'^A{x) 



401408 2 56 

1 H X X A(x), 

47403225 153 ^ ^ 



since C{x) > from ( p.24| ) and 9A{x)C{x) - AB{xY > turns out to be 



a(A) < 



1032757r 
896 



3.62 X 10' 



Positivity of f"(ro(a;); x) 

v''ir;x) = 2 



15C(a;)r^ - 6B{xy + A{x) + 



> 



(3.38) 



|^a;M(x)2 + —x^Aixy, (3.39) 
70227 ^ ^ 102 ^ ^ ' V ^ 



(3.40) 



(3.41) 



(3.42) 
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Figure 2: Squared dynamical mass of fermion m? / shown as a function of the four- 
fermion couphng constant (yf^A^/47r^ = G for several fixed gauge coupling constants 
a (A) = e^(A)/47r, obtained from Eq. (|3.44|) . 



is also guaranteed, since the sufficient condition 9B[xY — 15y4(x)C(x) < is fulfilled. 
[This leads to a similar condition as Eq. ( p.41| ).] 
The final form of the partition function is 



Z[x\ = Yl 2^°7r3exp[-t;(x)], 
{site} 



(3.43) 



with 



v{x) = Vo{x) +Vp{ro{x);x) + -\n[vl{ro{x);x)] 



1 2 



1 



1, 

+ - In 
2 



12C(x){3 



2A{x)ro{xy + 2B{x)ro{x)^} 



1 - A{x)roixy + 25ro(x)4 



(3.44) 



where Eq. (|3.37|) has been used. Differentiating this with respect to x gives us the gap 
equation, whose numerical solution is shown in Fig. ^. In comparison with Fig. |l|, there 
is a large deviation in the region where both x = m^/A^ and G = g'^h? /Ati"^ are small, 
and the mass is increased by higher order corrections. 
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For sufficiently small x, xA{x) ^ 1 or x ^ «(A)/47r, Eq. ( |3.44| ) can be expanded: 

,,1 1 525x4 8 / 7 2/3 

vix) = fn x H m X ' 

^ ^ "^^ 3 2 544 3 V255 J 

25 / 7 V^^ / 3n 1792 



+ — — x^/M 4tt- ^^^-Inx (3.45) 
3 V255y V"(^) 172125 ' ^ ^ 



448 2 / Stt 224 



H x^ 



1377 V«(^) 20655 
yielding the gap equation 



Inx^ +0(x^/3) 



''^^'^ - i-e- + xE,(x)- ^ ^^V'^ 4 



2x V255y 



4: dx G ^ ' 2x \255j 9x^/3 

25/7 1/3 f Svr 523543 , \ 

H — — Inx (3.46) 

9 V255y V«(A) 688500 7 ^ ^ 

224 / 37r 21103 , \ ^, 5/3, 

H X — Inx + 0(x^/^). 

1377 V"(A) 41310 / 



The stability condition 



1 d^v(x) „ _ 1 / 7 \ 



Elix) + + 



4 9x2 ; V255/ 27x^/3 

7 25 / 3n 2589043 



+ V255J 27x2/3 (,a(A) 688500 ^"""^ ' ^^'^'^^ 



224 / 37r 62413 

+ 



1377 \a{A) 41310 



Inx^ +0(x2/^) > 0, 



is fulffiled for X < 1 and a{A) < 2589043 ^- Therefore an a-independent solution exists: 

^' - 2(1^ <='-^*^> 

for X obeying x <^ Q;(A)/47r and Eq. ( |3.29| ). In this case g'^ is also highly suppressed: 

ttG = ^ 2x ^ 2. X 10"^^ (3.49) 

47r 

for X ^ 10-^2^ ^i^at is, A ~ 1 TeV. 



4 Discussion 

We give a gauge invariant recipe for calculating the effective action in QED with the 
four-fermion interaction. We use perturbation and the local potential approximation 
to study the dynamical fermion mass. In order to include the higher orders of the 
gauge coupling constant we just expand the completely known function ( p.5| ) into a 



15 



polynomial of + Fl[= F^^F^^/2\ and FlF^[= [F^^F^^jAf] and integrate with 
respect to the proper-time r. The last task is to evaluate the triple integral of s, 
t, and u. In this article we employ the WKB approximation for it, which, however, 
cannot be performed in an elementary manner when higher orders are included further. 
Meanwhile the ordinary perturbative treatment, considering the Gaussian part the 
kernel and expanding the exponential of higher parts than the third order, can always 
be ensured. 

As for the local potential approximation its efficacy is still open. In order to go 
beyond the local potential approximation, we must perform the functional integration 
of S', T, and [/, instead of the ordinary integration like ( |3.18| ). For example, an 0(a(A)) 
quantity after integrating with respect to the gauge fields would be 

(fx^ (fy^ {J^A^{x)JpAp{y)) / / (fx^ 



where 



l^^ip") {P\P-P.PP) J^,{FMF,A-P)), (4.1) 



1 

{Fp^{p)Fp^{-p)) = — {PpPpSy^ - PpPJup - PuPp^i^a + PuPJpp) , (4.2) 



so that 



/ - - 1) = /^n(p^) <^>-M^>--'-rt> . (4.3) 

The local potential approximation is to expand n(p^) in terms of p"^ and keep the 
lowest. Since the change of variables, Eq. ( p.l4| ), gives Fpi,{p)Fpy{—p) = S{p)S{—p) + 
T{p)T{—p) + U{p)U{—p), the leading correction to the local potential approximation 
is 



Jd^x n'(o)(-5n^5 - TD^T - un^u). 



(4.4) 



Therefore we must calculate the effective potential with the help of the Feynman graphs 
using the propagators of 5*, T, and U. 

Our result indicates that chiral symmetry is always broken for any a (A) and G as 
is seen from Figs. |l| and ^. This seems to be different from previous results H, |^, 
which claim that there is a nonzero Gc if a (A) < vr/S. It is, however, too early to 
conclude, since our present work is restricted only to the tree level of auxiliary fields 
and thus no higher orders of G are included. As is seen from the Figs. || and |^, the 
small negative correction, linear in G, to dynamical mass could easily swallow the 
broken region. From a recent study[|l^, the one-loop inclusion of the auxiliary fields is 
promising. Thus a further study must be necessary for a definite conclusion. 
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In QED the functional integral in terms of the field strength can be constructed, 
since the configuration space of the gauge potential as well as the field strength is trivial 
enough for the gauge to be completely fixed. It is challenging to generalize Eq. ( |2.33|) 
to QCD where an obstacle for gauge fixing, the Gribov ambiguity|jl8|, exists. In order 



to examine the dynamics of chiral symmetry breaking and color confinement in QCD, 
this must be done and the choice of convenient variables describing the low energy 
phenomena |jl9[| must be necessary. 

These directions of study are in progress. 
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